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Abstract 

A duality formula, of the Hardy and Littlewood type for multidi- 
mensional Gaussian sums, is proved in order to estimate the asymp- 
totic long time behavior of distribution of Birkhoff sums S n of a se- 
quence generated by a skew product dynamical system on the T 2 torus, 
with zero Lyapounov exponents. The sequence, taking the values ±1, 
is pairwise independent (but not independent) ergodic sequence with 
infinite range dependence. The model corresponds to the motion of a 
particle on an infinite cylinder, hopping backward and forward along 
its axis, with a transversal acceleration parameter a. We show that 
when the parameter a/ir is rational then all the moments of the nor- 
malized sums E ((S n /y/n) ), but the second, are unbounded with re- 
spect to n, while for irrational a/ir, with bounded continuous fraction 
representation, all these moments are finite and bounded with respect 
to n. 
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1 Introduction 

Diffusion processes are the most important transport processes to be mod- 
elled in the frame of dynamical systems. For a particle having deterministic 
"chaotic" motion, diffusion appears in the long-time limit of the displace- 
ment governed by the dynamics under invariant probability measure. The 
extension of the Central Limit Theorem to deterministic dynamical systems 
is well established for many strongly chaotic systems. In these systems dif- 
fusion is related to strong mixing properties of the motion of the particles. 
However, this condition is sufficient, but not necessary. 

A stationary bounded sequence of identically distributed random vari- 
ables Xq, Xi, . . . , with zero mean value is ergodic if: 

1 n— 1 

-Vl, — > E(X ) = (1.0.1) 

n *■ — * n— »oo 

i=0 

It fulfills the Central Limit Theorem if there exists a positive sequence f n 
such that j- Y^=o converges in law to a normal centered distribution with 
variance 1. The validity of the Central Limit Theorem goes beyond the case 
of independent (Bernoulli) sequences Xj. Recently this problem has been 
cast in the more general setting of ergodic dynamical systems. 

Let T be a measurable transformation on a space X, equipped with a 
probability measure pna cr-algebra T of measurable sets. The measure \i 
is supposed to be T-invariant: for all A of T one has /i(T~M) = ^i(A). The 
system is also assumed to be ergodic: 

n-i r 
V/e^W; -J2f(T n X)^ f(x)d^x) = E(f) (1.0.2) 

almost everywhere. The sum: 

n-1 



S n f(X) = J2f(T n X) 



i=0 



is called a Birkhoff sum. A real square-integrable function /, with zero 
mean value, is said to satisfy the Central Limit Theorem if: 

lim V [X : -== £ f(T n X) €/=-=/ exp(-^)rfx (1.0.3) 
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In other words, the sum 



y/«S n f>} 



= Y17=o f°T n converges in law to a Gaussian 



variable of zero mean value and variance 1. A generalization of the Central 
Limit Theorem has been established for a class of regular functions for K- 
systems [Go], and for hyperbolic dynamical systems the auto-correlations of 
which decrease exponentially (see a review in [Li]). One of the most remark- 
able results of the theory of chaotic dynamical systems was the convergence 
of normalized Birkhoff sums of observables to a diffusion process. In the 
Lorentz gas, Bunimovich, Chernov, Sinai [BCS] , [BS] and Young [CYo] have 
proved such behavior. They found a class of functions / for wich the rescaled 
ergodic sums converge to the Wiener process, i.e. 



where t G [0, 1] and Y t is a Brownian motion (other informations can be 
found in [Bi]). In non integrable area preserving maps numerical results on 
diffusion have been obtained by many authors (see, for instance, references 
in [Za], [MMP]). However to our knowledge, in the hamiltonian dynamical 
systems no rigourous proofs have been given. In [Za], there are several nu- 
merical evidences of the existence of anomalous transport in hamiltonian non 
integrable systems on account of the existence of islands into islands in the 
phase space of the systems due to the abundance of some stable periodic 
behaviors generating Levy flights. One of the main motivation of our work is 
to investigate diffusive behavior in simple ergodic area-preserving mappings 
with weakly chaotic properties (i.e. randomness with zero Lyapounov ex- 
ponenets [CH3]). In such systems, the problem is to find a class of functions 
/ for wich the rescaled ergodic sums converge to the Wiener process. A first 
step toward such result is to estimate the asymptotic behavior of the Birkhoff 
sums. In this work we consider this problem for completely non-hyperbolic 
systems having nevertheless a family of observables for which the two-times 
auto-correlations are zero. Swante Janson [J] has studied examples of non- 
ergodic stationary sequences with zero mean value and variance a, such that 
the sequence S n — Sj=o -^-i converges in law, without normalization, which 
implies a breaking of the Central Limit Theorem. 

We shall consider an ergodic sequence of pairwise independent variables, 
each of them taking values in { — 1,1} and having variance a 2 = 1. They 
are derived from a dynamical system on the torus T 2 , called the Anzai skew 



[t/r] 




(1.0.4) 



n=0 
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product and defined by: 



T(x, y) = (x + a,y + x) mod 2n 



(1.0.5) 



When a/n is irrational, Furstenberg (see ref. in [CFS]) has shown that 
this transformation has a unique invariant measure, which is thus ergodic 
and is in fact the normalized Lebesgue measure In a previous work 

[CHI] and [CH2] have determined partitions {^4,A C }, which are pairwise 
independent: n(AC\T~ n A) = (J>(A) 2 . It is equivalent to say that for the func- 
tion f — I A — fi(A) the family of functions / o T n are pairwise independent 
random variables. However, this does not imply at all that these variables 
are jointly independant. Actually, the process has infinite memory because 
its metric entropy is zero, moreover, there is only linear divergence of trajec- 
tories. It is to be noted that on account of the pairwise independence, it can 
be immediately seen that the "diffusion coefficient" a defined by: 



is finite. However, the existence of the diffusion coefficient is not sufficient 
to imply asymptotic normal distribution. 

Thus it would be interesting to see whether the Central Limit Theorem 
holds here or not. We shall see that the result depends essentially on the 
continued fraction representation of a. In this work, A will be limited to the 
subset Tx [0, 7r], Xi = f oT l and the moments E((S£)) will be asymptotically 
estimated. It is however difficult to reconstruct the characteristic function 
$ n (t) of ^js=. For odd q, E(S£) = 0. For even q several situations can occur: 
We shall show that: 

a) If a/ir is an irrational number with a bounded continued fraction rep- 
resentation, then E((S 2q ) = 0(n q ). Consequently | E((^) 2q ) \< A q , where 
A q is a constant. The convergence of the characteristic function depends on 
the A q . 

b) If a/n is an irrational but not of the previous type, then there is no 
/„ such that E((j^) 2k ) = 0(1), different from zero for all k > 0. 

It is to be noted that our mapping gives interesting examples where the 
existence of th diffusion coefficient a 2 = lim^^ E(S 2 )/n = 1 is not sufficient 
to imply the convergence to the normal distribution. 




rt-l 
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So it may be suggestive to start out with a concrete picture of particle 
dynamics. The dynamical system Eq. 11.0.51 corresponds to the motion of 
a rotator kicked at regular time interval by a force modulated so that the 
angular velocity remains bounded. Let 9{t) and 9{t) be the angle and the 
angular velocity of the rotator, where 9{t) G [0, 2n[. Then, at time t + 1, the 
state of the rotator is given by 

9{t + 1) = 9{t) + 9(t) mod (2tt) 

0(t+l) = 6(t) + a mod (2tt) (1.0.6) 

and we have at time t: 

9{t) = 9(0) + t6(0) + ^~ 1 ) a mo d (2tt) (1.0.7) 

An idea of the randomness of this motion can be seen by using a partition 
{A,A C } of the torus into two regions: A = {9 e [0,7r[} and A c = {9 6 
[7r,27r[}. It can be seen that: 

fi(Ai n T~ n Aj) = n(Ai)n(Aj) 

for any n ^ 0, where A t is either A or A c [CHI]. This partition is far from 
being a Bernoulli one since the entropy of the system is zero. We call such 
partitions pairwise independent. In [CH3], it is shown that such sequences 
are unpredictable in the sense of Wiener least squares criterion. It is natu- 
ral to study the distributions of sums of such sequences, which represent a 
particle displacement induced by the dynamical system of the Eq. (jl.0.5|) as 
follows. 

A particle moves on an infinite plane among periodically distributed ob- 
stacles with spatial period equal to 1 along both (qi, g 2 ) -directions. In the 
gi-direction, the motion of the particle is uniformly accelerated at each reg- 
ular time interval by an amount a and has uniform free motion along the 
(^-direction. That is, define the velocity pi(n) — qi(n + 1) — qi(n) , then the 
equations of the projection of the motion in the gi-direction are : 

qx{n) =g x (n-l)+pi(n-l) (1.0.8) 
p x (n) =pi(n - 1) +a (1.0.9) 

It is the result of the discrete time action of the mapping: T : (pi,qi) — > 
T(p 1 ,q 1 ) given by: 

T(p 1 ,q 1 ) = ( Pl + a,q 1 +p 1 ) (1.0.10) 
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Figure 1: Trajectory of the particle in the plane ((71,(72)- 
Thus, we obtain: 

qi (n) = gi(0) + npi(O) + n(n - l)a/2 (1.0.11) 

The particle is moreover submitted at the begining of each time interval to 
a field changing the direction of the motion up and down along g 2 -direction 
in terms of its position along the gi-direction in the following way (see figure 
1): the velocity direction of the particle p 2 at time t = n is given by x(<?i( n )) 
where x( x ) * s a periodic discontinuous function defined by: 

xW = /-l if ^]0 1/2] 

AV 1 \ 1 if x e]l/2, 1] v y 

That is, the direction of the velocity at time t — n is : 
e„(pi(0),gi(0)) = x(qi(n)) = X (?i(0) + npi(0) + n(n - l)a/2) (1.0.13) 
and the value of the variable g 2 at time t — n + 1 is: 

n 

g 2 (n + l) = g2(n)+x(9i(n)) = 5]e<(Pi(0)>9i(0)) (1-0-14) 
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Let us denote by Sk the position of the particle along the (^-direction at 
t — k. The equation for such a motion is expressed by the following recursion 
relation: 

S n +i = S n + x(qi{n)). 

S = (1.0.15) 

In terms of the variables (x, y) G T 2 the distance travelled by the particle 
along the vertical direction after n steps is: 

S n (x, y,a) = J2 X(y(k)) = J2x(y + kx + o\ (1.0.16) 

The problem is to study the limiting distribution of such random walk. 

In the following sections, the proof will be given in steps with an increas- 
ing order of complexity, which turns out to be related to number theoretical 
questions. The main objective is to obtain the asymptotic behavior of the 
expectation values of S^(a) for n — > oo. In the course of this study, we 
encounter the so-called multidimensional Gaussian sums which are a par- 
ticular case of the Weyl sums [K]. They are essentially sums of the form 
used in the definition of theta functions, except that the summation is finite. 
The Gaussian sums fulfill an exact duality formula, known as Landsberg- 
Schaar formula [L] , which is a variant of the famous Jacobi identity for theta 
functions. Except for the rational assumptions in the application of Schaar 
formula, some works had been done to estimate the growth rate of the one 
dimensional Gaussian sums by Hardy and Littlewood in [HL]. To this end 
they have first established a duality formula and then used the theory of con- 
tinuous fraction decomposition (here we use tools and notations from [VdP] 
and [R]). We shall treat respectively in sections 3, 4 and 5 Gaussian sums in 
one, two and d dimensions, the central point in each part is concentrated in 
the proof of a duality formula for these sums. 

The Gaussian sums have been more recently studied from a geometrical 
point of view. [DM] and [Z] have applied the results of [HL] on polynomials 
associated to modular functions. Interestingly one may derive this identity 
from a simple quantum mechanical system having the torus as phase space 
[AR]. But Gaussian sums are found also in other fields like the problem of 
fractional wave packet revival [SLB]. 
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2 Preliminary considerations and the case of 
a rational 

Before going into the crux of the subject, let us introduce some technical 
steps. The important quantity to study here is the q th moment defined by: 

E(S? l (a))=[ (S*(x,y;a))^ (2.0.1) 

An objectve is to study the convergence of the characteristic function: 

$ n (f) = f e it.(S n ( X ,y-,a)/V*)^y = y M! jB (^( a ) /v ^) (2.0.2) 

JloMxioM A1X ^ 9 ! 

Using the Fourier decomposition of the function x' 

n p iPy 



where Z is the set of relative odd integers, we have an alternative form of the 
distance travelled by the particle: 

s n (x, y;«) = ^EE p • ( 2 - - 4 ) 

Hence with this form one may compute the q th moment. After integration 
on x and y, which yields two constraints: 

X> = o 

i=\ 

Q 

Y^hPi = (2.0.5) 

i=i 

The above expression becomes: 

E[S<{a)] = (f V E • • • E E • • • E ^-pe^^-^ 

(2.0.6) 

where (5 — We are now in a position to prove: 
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Lemma 2.1 

E[S^\a)] = 

M.(-^£...££...£ ^i 

x exp J2 E v^V^ - ( 2 - - 7 ) 

j=2 j=2 <^Z=2 P 

where 5(.) is the kronecker function: equal to 1 when the argument is and 
otherwise. 

Proof. (Pj)i=i,...,2g are odd. The condition X^i=i Pi — ^ can only be satis- 
fied if q is even. Hence E[S^ q+1 (a)] = 0. Now from the constraint X)i=i = 

0, we get Pi = - X]^ 2 p i and from Ei=i ^Pj = i- e - Ya=2 Q P i( k i ~ ki) = °- 
Thus, working out k±: 

J2i=2 Pjk 

E 2 i p 
i=2 

with (&i)i=2,...,n-i ^ {0, ... ,n — 1}. And we obtain, by substitution 



fci = ^= 2 (2.0.8) 



2g / 2g 2g 2 g 2g 

i=l 2^i=2 \ j=2 j=2 i=2 j=2 

Z^i=2- r « V i=2 j=2 / 

which ends the proof. ■ 

In particular E[S^(a)} can be exactly and directly computed : 

p i£ z p 2 ez fci=o fc2=o 

But we get from the constraints, Pi = — P 2 = P and ki = k 2 = k. One has 
thus: 

= ^£^£i (2.0.H) 
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From the property of the Riemann Zeta-function [MOS], it may be shown 
that Y.p& l / p2 = ^V 4 - Hence 

E[S 2 n (a)]=n (2.0.12) 

The fact that odd moments are zero implies that the characteristic func- 
tion is also even and consequently the Fourier inverse transformation of the 
characteristic function, that is the probability distribution function of the 
random variable S n (x,y; a) is also even. 

Finally we observe that the series 

£ (2013) 

(P 2 ,..,P 29 )ez%-i 2gK 2 ^ 2q> 



is absolutely convergent since, for all (P 2 , . . . , P2 q ) £ Z 29_1 , \P 2 + . . . + P 2q \ > 
1. Let us prove it for q = 2. The generalization to arbitrary q is straighfor- 
ward. The series may be split into: 

y I = 2 y I 

(P 2 ^3 \^ P ^ + * + P 4 )l {p2 ^ 3 | + PS + P<)\ 

, ^ . \p 2 pzp±(p 2 +pz-p±)\ 

(2.0.14) 

Let (x, y, z) G [0, oo[ 3 and < 8P < 1 such that: 



1 < P 2 < x < P 2 + 5P x - 5P < P 2 < x 
1 < Pz<V < P-3 + 5P y - 5P < P 3 < y 
1<Pa<z<P a + 5P <£> z-5P <P A <z (2.0.15) 

The first sum in the r.h.s. of the equation 12. 0.141 is bounded : 
^ \P 2 PMP 2 + P 3 + PA\ ~ 6^3 ^ 

{P2,P3,P i )&i S (P 2 ,P 3 ,-P4)GN 3 
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The second sum can be transformed using: < (x — 5P) (y — 5P) (z — 5P) < 
P2P3P4 < xyz and x + y-z- 25P < P 2 + P 3 - P 4 < x + y - z + 5P. Hence 

< ——. — 1 —— — (2.0.17) 



|P 2 P 3 P4(P 2 + P 3 - Pa)\ ~ (x- 5P)(y - 5P)(z - 5P) 

1 1 



max 

def 



\x + y — z — 25P\ ' \x + y — z + 5P\ 

F 5P (x,y,z) 

But V5P < 1, 3e > 0, such that any hypercube 

[(P 2 , P 3 , P 4 ); . . . ; (P 2 + 5P, P 3 + 5P, P 4 + 5P)\ 

does not lie in the intersection il t d = [0, oo[ 3 n{|x + y — z\ > e}. Then we 
have 

^ ^P 3 P,(P2 + P 3 -P,)\ (2 - - 18) 



1 

< 



(5P) 3 

(P2,P 3 ,P4)eN 3 



rP 2 +5P rP 3 +8P rP 4 +5P 

E / / / Fsp(x,y,z) 

'p 2 Jp 3 Jp 4 



- 7F777T / F sp(x,y,z) 
'[o,oo[3\n e 



(5P) 



One can easily estimate the left-hand-side by realizing that the integration 
of the integrals: 



f dx.dy.dz f 

ifo,oof3\n f \x + y - z + 5P\ L 



dx.dy.dz 

(2.0.19) 



'[o,oop\n e \x + y - z + SP\ ' 7 [0 ,oo[3\n £ \x + y-z- 2 X P\ 
yields an absolutely convergent result. 

2.1 a = 

We consider the simple case a = and show that the case with a rational 
can be brought back to the case a = 0. First consider again: 

S n (x,y;0) = -J2^ elP{kX+V) (2- 1 - 1 ) 
in — i r — ' 

Pel k=0 
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The above geometric sum can be reexpressed in terms of u — e lx and v = e w 
as: 

P<=7, ^ ' e=±l 



Now, the 2k moment 




(2.1.3) 



may be computed as complex integrals evaluated on the the product of the 
circles of unit radius and centered at 0, i.e. C(0, 1) x C(0, 1). In the domain 
bounded by these circles, the integrand is biholomorph and the product of 
the contours is homotopical to C(0, 1/2) x C(0, 1/2). Thus (1 — uP) -1 can be 
expanded in power series around 0, with u = e tx et v = e iy : 



«.(«. »; o) = | E« f ( E I (E «"' p I ( E ^ v 

Pe Z \r?=±l / \m=0 / \e=±l 

(2.1.4) 

or 

S B («, v; 0) = | £ ^ ^u p (^^)v^ (2.1.5) 
We can now compute E[S^(0)] with this expression: 

e[^(o)] = -± y y mmi^**** 

77^,772— ±1 mj_ ,m2 €iV 
ci,e 2 = ±l P 1) P 2 eN 

(fa ^ Pl (n+zn n+mi+ i^i-) + p 2 (a+ZZ2„ +m2+ i^2-) 

2i?r 

(2.1.6) 

The first integral over t> yields eiPi + e 2 P2 = 0, one must have e-y = — e 2 and 
Pi = P 2 = P. The exponent of z in the 2 nd integral becomes: 
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A I — g — n + mi H 2 — j V — 2 — n + m2 H 2 — 

p ,'Vl+V2 , . 1 

P — - — n + mi + m 2 + 1 



(2.1.7) 

The contribution to the complex integration comes only from the vanish- 
ing of this quantity, this leads to the following: 

Vi = V2 = -1 (2.1.8) 
m i + m 2 — n—1 (2.1.9) 

and card{(mx, m 2 ) G IN 2 tel que m 1 + m 2 = n — 1} = n, thus: 

e=±1 peN Pez 

This computational method may be extended to higher moments E[S 2k (0)}. 
We see that using Eq 12. 1.61 to power 2k yields 

Eis n m = <-i)» (±V E E E E 

»7i,-»72fc=±lei,-«2fc=±l mil ... m2fce INi J i,...P2*6N 

2fc 

Al ^ / 2i7r / 2ivr v • • ; 

We can perform the t>- integration and get: 

£^ = (2.1.12) 

4=1 

and the ^-integration yields the following proposition: 
Proposition 1 

£[s„ ( o)"] = (-d'(^)' E E E E 

r ? i r ..7?2 fe =±le 1 ,...e 2 fc=±l mir ..m 2k eIN Pi,...P 2k eN 

2k 

W^M^e^Pm) (2.1.13) 

i=i 1 
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where the coefficients A n are defined as: 

2k , 2k 



A n [ri u e h Pi](2k) = card{(rrn . . . , m 2k ) e W 2k | ^ rmPi = -- P(l + nrji)} 

i=i i=i 

(2.1.14) 

w^/i £/ie constraint Yliti e %P^ — 0- 

This result is obtained by the Residue Theorem. It remains to determine 
the behavior of A„ as n — > oo. We observe that not all the combinations 
of e,- are to be taken into account since (Pi, ... , P„) G N 2fc . Thus the choice 
(ei, . . . , €2k) = (1, • • • , 1) is excluded by the constraint. According to the 
usual notation , if the (Pi)i=i,..., r do not have common divisors, we note: 
(Pi,...,P r ) = 1, then: 

r i n r ~ l 
card{(m 1 ,...,m r )eW\y^m i Pi = n} ~ — — — (2.1.15) 

■J n-+oo [r — P± . . . P r 

Before going on we prove the following lemma: 

Lemma 2 Let YZ,=i = an< ^ (-^ >1 ' • • • ' = 1 * s equivalent to (Pi, . . . , 
Pj, . . . , P r ) = 1 where Pj, is removed. 

Proof. 

It is clear that if (Pi, . . . , Pj, . . . , P r ) = 1 then (Pi,P r ) = 1. Actually if 
(Pi, . . . , Pj, . . . , P T ) = 1 then fls E N \ {1} such that Vi ^ j, P = sp* and 
thus (Pi,...,P ? ) = 1. 

Conversely if (Pi, . . . , P q ) = 1, but (Pi, . . . ,Pj, . . . , P r ) = s with s e 
N \ {1}, i.e. 7^ j, P = spi, then Pj = Y : .. j< ,<;,!'; = ^jPi) whi ch 

contradicts the fact that (Pi, . . . , P q ) — 1. ■ 

Coming back to our estimate of A n , put IN 9 = N 9 n {(Pi, . . . , P q ) = 1}. If 
(Pi, ... , P2k) e IN 2 " such that E?=i uPi = then: 



i=i 



1 if ^2i=iViPi — ® (2.1.16) 
otherwise 
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with 5^i=i e «-f« = 0- Now from Eq l2.1.lfll the asymptotic behavior for n — > oo 
of the 2/c th moment is: 



B|s »' (0)l ™ 2 E ^ E e n f ( E 



P £ N <«=±i p ^e# 2fc i=l * \i=ia 
H =±i 

(2.1.17) 

with 52i=i e «-^ = anc ^ Si=i ^i-fi > • We have already seen that the series 



2k /2k x 2k ~ l 



E E n^(E^.) 

f7i,...,»72fc=±l P 1 ,...,P 2k ^N 2k i=1 1 Vi=0 / 

is an absolutely convergent series. Since E^opS > 1, we have thus the 

(J— >oo 

theorem: 

Theorem 2.2 VA; G ZZV, (0)] ~ C k n 2k ~ x with 

d = l 



^^Enf e^< ("■») 

v y li.tpii i=l 1 \i=l / 

2.2 Extension to rational /? ( a £ Q) 

To this end, we regroup terms of S n (x, y; a). Wa G Q, 3k G IN such that 
VA; < k : 

k(k-l)P ^ mod (1) 
A; (A; - l)(3 = mod (1) (2.2.1) 

We shall set n = qko and write S n (x, y; a) grouping the terms such that: 

3-1 

5 s (x, y) = J2x(y+[s + nk ]x + 2s(s - 1)tt/3) (2.2.2) 

n=0 



2 THE RATIONAL CASE 



16 



and 

fco-l 

S n (x,y;a) = ^26 a (x,y) (2.2.3) 

with < s < k . The sum can be partially performed after inserting the 
definition of x- 

S s ( x ,y) = — > — e * fe ° 2 x — ^ 

Pez 

^p(j/ +a; [s+M|zi)] +2s(s _i )/3) _ e - iP{y+x[s+ M*pV ]+2s{s _ m ^ 

Call now s' = 2s(s - 1), u = e", u = e i?/ and t = e 2 ^. Then 
5 s (x,y;a) = 5 8 (u,v;t) = - ^ -u * 

ez tP(s + ^k ) v eP t ePs' (2>2>5) 

c=±l 

We now expand j^p^ as in the previous section, and obtain consequently: 



6s ^ V]t) = ± J2 ^ fP f Ps u p ( fS+fc » 1 ? +mfc » + * 2 f) (2.2.6) 

e,r,=±l PG 2, me ]N 



The expectation value: 

' %,e,;=±l 1=1 * 



PiSN.m^elN 



(2.2.7) 

is not zero if Y^Li Pi{ e i s i + ^o^r 1 + mko + g&o 2 ^) = and Ynti e i p i = °- 
As before, one can reduce to the case (Pi, . . . , P2k) = 1- The first assertion 
is true if Yli=o €iP>iSi * s divisible by k , with < Sj < fco- The hypothesis on 
Pj suggests that there exists a P io not divisible by fc . 
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Lemma 3 Let s±, . . . , s 2 k such that X^=o e «-^ s « ^ s divisible by k then Vs e 
{1, . . . , k } and s ^ s io , k Q is not a divisor of Y^t=o e iPi s i- 

Proof. Suppose that this assertion is not true: 3n such that Y^=o e iPi s i = 
nko and 3n' such that Y^=o e iPi s i = n 'ko- Hence ej Pj (sj — s) = (n — n f )ko 
but k is not a divisor of Pj and |sj — s\ < k , this is not possible. ■ 

Consequently card{(si, . . . , 82k) G {1, . . . , ko} 2k } such that ko is not a divisor 
of Ylfti e iPi s i an d the value (Pi)i=i,...,2fc are mutually prime numbers of values 
less than k^' 1 . 

Remark 4 

We have 



card 



2k 



i=l 



(mi,...,m 2fc ) G IN 2fc I ^P(e^ + /cq 1 6t + m^ + qk ~ ^ ^ ) 



2fc 2fc 

9 



~ card[(m 1 ,...,m 2fe ) G 1N 2A; |V m,P = -J V^P] (2.2.8) 

i=l i=l 

with (P 1; . . . , P 2k ) = 1. Since ^ = and t p ^=^ PiS 'i = t p ^ eiPiS i/ 2 . 
Following the previous procedure we obtain for n = qk . 



2k 



E[S$ (a)] ~ 

g^oo 



2k 



,i=l 



v 7 ^=±1 i=i 



2fc-l 



1 



k. 



2k-l 




pgN 



Si = l 



p2fc 



(2.2.9) 



with the conditions Ylfti e i^i = >zCi=i e i-P« s i is divisible by /c , (Pi, • • • , P2q) 
1 and 5^i=i ^j-Pj ^ 0. As it was mentionned , for all P G N, we have 



< 1 



(2.2.10) 



Note that for a = (fco = 1) this term is simply I. Consequently we have 
the theorem: 
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Theorem 2.3 Va G Q, \/k G IN* E[S* k (?f)) ~ C k n 2k ~ x with d = 1 and 
Ck a constante. 

This result suggests a unique normalization, for which the characteristic 
function of S n (x,y; a) once normalized will be independent of n, when n — > 
oo: 

Consequently the probability distribution function F of the normalized ran- 
dom variable limit S = ]ini n ^ 00 ^Mi£^ ; for a G Q: 

F{S) = S (2.2.12) 

where 5 is the Dirac distribution. That follows in fact from the ergodic 
Birkhoff theorem, but the speed of the convergence to zero of i-MEilig) j s very 
slow. 



3 The case of irrational (3 in terms of Gaus- 
sian sums 

In this section we shall proceed with a different type of estimate of E(S^ q (a)) 
for a G IR*. In particular, we emphasize the difference between the rational 
case of section 2 and the irrational case a G IR \ Q. For this purpose, we 
shall establish a generalization of the duality formula for one-dimensional 
Gaussian sums, given in the celebrated work of Hardy and Littlewood [HL] . 
This problem consists in estimating the behavior of the many variable sums, 
for large n: 

n-1 n-1 . asr 2 q V 2„ PjPj n , v 2 

e ?EwE,rf ^ ^ (3.0.1) 

fci=0 k 2q =Q 

where Pi G Z, i — 2, . . . , 2q are parameters (see Eq J2.().7j) . The above sum 
bears over integer in indices of IN 29 formed by intersection of the hypercube 
of total volume n 2q and the hyperplane defined by the equation: 

2q 

J2 p i(ki~h)=0 (3.0.2) 

i=2 
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This condition can be recast as an integral using the formulae: 



i 

2-nixP , 







e Zmx ^dx = 1 if P = 

= if P G IN* (3.0.3) 

Define now, the quantity: 

£ n [x,/3;P 2 ,...,P 2? ] = (3.0.4) 

fci=0 fc 29 =0 

Then, the 2g th (see Eq J2.().7j) moment may be rewritten as: 

E[S%A = (-I)* (-Y (3.0.5) 



4 V 71 ", 

1 Z" 1 

X fr2g p /^2g p , / A P 2 , • • • , P 2g ]ofe 

As £ n [x, /3; P2, . . . , P2 9 ] are analogous to 2g variables theta functions [M], 
we shall express the formula Eq. 13.0.51 in terms of Gaussian sums d crj^ 6 '(f2, 9). 

Definition 5 we define the d- dimensional Gaussian sums: 

d af^(n,9) = X e^-l^-^^H^IW-l) (3 . . 6 ) 

Where Q is a d square matrix with real coefficients: Q G Md(IR), 9 G 
Vectd{IR) such that Vi = l,...,d, the component 9i of 9 is restricted to 
< 6i < 1. Here we denote by t v the transpose of a vector v. Moreover 
a and b are elements of Vectd({0, 1}). Af t , with n = (m, . . . ,nd) , is the 
hyper-rectangle with integer sites in IN* defined by: 

Any k G A^ is such that \/i = 1, . . . ,d, < ki < rij — 1, ki G M. 
With this notation: 

n-l 

e n [*,/3;P 2) ...,P 2? ] = ^ e - 2! ™E^^^-V?(n,0) (3.0.7) 

fci=0 

where: 
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• An defined by the 2q — 1-vector ft: n — (n — 1, . . . , n — 1) 

• the 2q — 1-square matrix Q: 



n 



P 



P2(J2-i 2 Pi-P2) -P2P3 

-P2P3 PsiT&Pi-P*) 



-P 2 P 2q 



P 3 P2 q 



• the 2q — 1-vector 9: 



9 = x 



P2 
P 3 

P 2q 



E Vect 2q -i(IR) 



— P 2 P 2q 
-P 3 P2 q 



P2cS52?£L 2 Pi - P2 q ) 



(3.0.9) 



• a = b = 

In order to compute the above sum, we shall first establish the duality for- 
mula: 



(3. 



Theorem 3.1 (duality formula) IfQ e M d (]0, l[) sym ,inv ( the d- dimensional 
symmetric and invertible matrices with values in ]0,1[), 6 e VectdQO, 1[), and 
(a,b) e Vect d ({0, 1}), then 



d <j$ b] (n,9) = 



ri d d [b,S] . 
• Vi 1 ■■■Vi d (CIV) 



(3.0.10) 



Vdetfi/ 

• ' \0<ii,...,t c i<a 

where T> is a domain of IN d and Q(T>) is the set of integers sites of IN d 
contained in the image of V by Q and V{ is the projection on {k; L = 0}. 
Moreover we take the convention 0(1) = 1. 

This duality formula will be extended to fl G M d (JR) sym;inv (d-dimensional 
real symmetric and invertible matrices) and later we shall treat the case of 
singular f2, i.e.: det(fi) = 0. 
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3.1 Results of Hardy and Littlewood on one-dimensional 
Gaussian sums 

Following [HL], let us introduce: 



n-l 



(%(x,6) = J2e inx(k ' 1/2)2 cos(2k-l)n9 

k=0 
n-l 

C*(x,9) = ^ cos 2k7ie 

k=0 
n-l 

C*{x,9) = ^(-l)V 7ncfc2 cos2/br# (3.1.1) 



k=0 



Then, the duality formula for the Gaussian sum C^(x, 9) (see Eq. 13. 1.71 below 
) can be obtained by application of the residue theorem in the formula [HL] : 



C 6 n {x, 9) = W - e -™-C 6 nx {— , -) + i / . r - Q»(fl, x, t)dt (3.1.2) 
V x xx J sinh7rt 

with < x < 1 and < 6 < 1. The function Q n (#, x, t) is given by: 

Q n (0, x, n) = cos(2vm#) cosh(27r0t) sinh((2nx -2k- l)nt) (3.1.3) 
+ i sm(2nn9) sinh(27T0£) cosh((2nx - 2k - l)nt) 

Observe that the integral in this equation may be reduced to integrals of the 
type : 

sinh(Trt) v y 

From the theorem on intermediate values we have: 3T bounded, such 
that V# g]0, 1[, the previous integral is 

pT(0) 

Jo 

where 0(1) means a bounded function with respect to 9. This integral de- 
pends on 9 but remains bounded. The validity of this theorem is guaranteed 
since ^fer is a positive decreasing function of 9 on IR + . Its value is: 



0(1) J -e~*£ ( exi{TVmt - - erf(-£=) ) (3.1.6) 

V x \ Jia Via 
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where the function erf (the error function) is defined as erf(z) i— > J* e~ t7TV dv 
This leads to a more precise form of the result of Hardy et Littelwood. A 
similar proof can be performed for C\ and C*. In fact we have obtained the 
proposition: 

Proposition 6 We have the following duality formulas: VO < x < 1 and 
VO < 9 < 1 

In their work, Hardy and Littlewood, have stated that analogous sums with 
sine instead cosine obey also duality formulas which can be established in 
a similar way. However as already seen, the parity of the cosine function 
does play an important role. When the cosine is replaced by the sine, there 
appears Fresnel functions which need additional treatment. We offer here an 
alternative proof. Define: 

q+r 

Cl q+r (x, 9) = J"** cos(2kn9) = C 3 q+r (x, 9) - C 3 q (x, 9) (3.1.8) 

k=q+l 

Applying the duality formula on C 3 (x, 9) we find: 

Cl+r^ 0) = yi e -*£ (cl qMq+r) (-1, ^ + 0(1)) . (3.1.9) 
But we also have: 

q+r q+r 

Cl q+r (x, 9) = cos(2ttA;#) = e lnx(k+q)2 cos 2n(k + q)9 

k=q+l k=l 

(3.1.10) 
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after expansion of the cosine, we obtain: 



-nxq 



2 r 



Cl Q+r (x,6) = ^ e ^ 2 (cos27T^ 



(3.1.11) 



k=l 



[cos 2n(xq + 9)k + cos 2n(xq — 9)k + i(sin 2n(xq + 9)k + sin 2ir(xq + 
+ sin 27r^g 

[sin 2n(xq + 6)k — sin 2i\{xq — 9) k — i(sin 2n(xq + 6)k + sin 2n(xq + 0) A;)]) 

Assume now x G Q, such that x = a/6, where a and b are prime numbers. 
We choose now e\ as multiple of b and let e(xq) label the integer part of xq: 

C 3 q , q+r (x,0) = (-1)**°*> (3.1.12) 

) — 1 r— 1 

cos 27rg# ^ e^ fe2 cos 2tt0A; - sin 2^ ^ e J7rafc2 sin 2tt0A: 



fe=0 



k=0 



Using the duality formula we obtain: 



CL +r (x,9) = (-iy^)cos2nq9J-e- i ^(C 3 xr (--,-)+0(l) 



x 

r-l 



X X 



.(_!)9e(*9) sin 27[q 9 e lwsk2 sin 2n9k. 



(3.1.13) 



fc=0 



Moreover we have 



x 



Ct 



---+o(i) 

x x 



1 p -™ e - 
— e x 

x 



e(xr)— 1 

E 

fc=0 



X 



(fe+e(^^)) 2 „ fc + e(xg) „ 

COS27T 



, . e(xr)-l 

cos 2% -———9 y: < 



fc=0 



fc 



, / \ e(xr)-l 2 

cos27r-0-sin27T^^0 V e- i7r Vsin27r-0 + O(l) 

fc=0 



(3.1.14) 
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This leads to the equality: 

e(xr)— 1 



COS 



= <£. ( -1 cos 2.^ - ^ f =i ^ sin 2.^ + (9(1) 

xr \ X X J X \ x xj X 

(3.1.15) 

when g is given as a multiple of b, we get simply: 



S r 3 (M) 



s„ 



-1 



x a; 



0(1) 



(3.1.16) 



M9 such that 2nq9 ^ irk or 9 ^ ^ mod (1). Let us now call 



C = jfl G]0, 1[ 3A; G IN such that = ^ mod (1) j (3.1.17) 

But V0, 3g' a multiple of 6 such that 9 (jL C?'- 

Hence, V# g]0, 1[ and Vrr G Q, we obtain the same duality formula . 
Finally S r (x,9) is a continuous function and bicontinuous since it is a finite 
sum of bicontinuous functions. As Q is dense 1R. Vx g]0, 1[, 3x n G Q\]0, 1[, 
n G IN such that lim^oo x n = x and by continuity we obtain the limits: 



lim 5 r 3 (a; n ,i 



lim S xr ( — , — 



1 



9' I — t 

X X 



(3.1.18) 



Doing the same proof for S^(x : 9) and S*(x,9) we get the following duality 
formulas: 



S*(x,9) 



-e'^ { St 
x 



x 



"r^r / C<3 

nx 



-,-) 0(1) 

X X 



1 



% 

X 



x [ S n x[ 6 



X X 



(9(1) 



(3.1.19) 
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These results are used to show that, for the following Gaussian sums: 



n-l 



1 a £' b \x,9) = ^ e -^( fc -«/ 2 ) 2 e 2 ^( fc - a / 2 )( e - b /2) (3.1.20) 



fc=0 



with a = 1,0 and b = 0, 1, we have Vx e]0, 1[ and 9 e]0, 1[: 

1 



Vj^(x,0) 



\ u nx 



X X 



0(1) 



(3.1.21) 



we summurise these results in the following: 
Theorem 3.2 Vx e]O,l[,V0 e]0, 1[, Va, 6 = 0, 1: 



Vj^(x,0) 

cM(x,e) 



-e VN 



<* 6] (M) = tk-Wt 1 (3.1.22) 



1 



x x 



1 



X X 



i e 



-0(i) 
0(i) 



X X 



Upon application of duality Vx e]0, 1[, V0 e]0, 1[, Va, 6 = 0, 1, we have: 

V^(-x,#) = 

-0(1) 



I 
x 



1 

x' x 



-e x 



'—x 



1 



Vj^M— ,-1+0(1) 



1 [6,o] 

^ — r) t 



X X 

1 

— x' — X 



+ 0(1) 



(3.1.23) 



we extend the validity of the duality formula to x e IR \ IN. Using the 
convention V - 1 = the identity formula can be extended for — 1 < x < 1 
and x 7^ 0. All this remain valid for C n and S n , since we have the identity: 



1 at b] (x,9)=Ct b \x,9)+ l .St b] (x,9) 



(3.1.24) 
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Finally we also have the formula: 

l a^ b \x + l,6) = Vi ai cr [ :> b ' ] (x,9) (3.1.25) 

with b' = b + a mod (1), it may be used to extend the validity domain of the 
formula into x G IR \ 7L. If x G 7L, the series reduces to a geometric series. 
The same results hold for C n and S n . 

Remark 7 

If x G Q, there exists a discrete duality formula (Shaar's formula). This 
relation is actually the same as our formulas without the term 0(1). The 
formulas are used in successive iterations until the summations are limited to 
a small number. The domain of summation is reduced at each step. However 
the result depends on its expansion as a continuous fraction (see article Hardy 
and Littlewood). This last problem is closely related to the Gauss transform 
T(x) = ^ mod (1) for which there exists infinite generating partitions. So 
for example if x is irrational, but admits a decomposition into continuous 
bounded fractions we have: 

1 a£> b \x,e)^o(VH) 

a homogenous function in 9. But if x G Q: 

1 a^ b \x,0)-O{n). (3.1.27) 



(3.1.26) 



4 Two-dimensional gaussians sums 

4.1 Two-dimensional Gaussian sums with non-degenerate 
matrix O 

Once this case treated, one can extend the results to Gaussian sums of ar- 
bitrary dimensions, provided some delicate points are worked. Consider the 
Gaussian sum in two dimensions: 

91-1 q 2 -l 

2 0^°\VL 6) = ^ ^ e i7r(wiA;2+w 2 fc2+2wfcifc2) e 2j7r(fci6»i+fc2e2) ^ ^ 

k 1= 0k 2 =0 
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where q = (q±, q 2 ) and Q G M(]0, l[) sym with 



UJi UJ 
UJ UJ 2 



(4.1.2) 



Af labels integer sites of rectangles of dimensions q\Xq 2 . We first assume 
that < ujq-y + 61 < 1 and < ujq 2 + 9 2 < 1. As q 1 and q 2 are to be large 
numbers, one may call these matrices weak coupling matrices. If uj — 0, this 
Gaussian sum goes back to a product of one-dimensional Gaussian sums. In 
general: 

91-1 92-1 
2 a^°\VL 6) = e i7r(a;iA: i +26,lfel) ^ e in(,u 2 kZ+2wk 1 k2+28 2 k 2 ) ^ j 3^ 

fel=0 k 2 =0 

We shall establish the duality formula successively with respect to (k 1 , k 2 ). 
If the duality transformation is performed on the k 2 sum one gets: 



fci=0 



!/ "'" J (o^fcx+fla) 2 




> e "2e "2 +0(1) 



CJ 2 



fc 2 =0 




a 2 



0)2(92-1) 

E 

fc 2 =0 



— e "2 1 } e "2 g w 2 e ^2^2 

C^2 



91-1 2 

»7r " ia,2 ~" fc 2 

g ^2 1 



E 

fci=0 

+ ©(^JXe-^^^-^^^e^-^)*! (4.1.4) 

* " fci=0 

when f2 is invertible matrix (det(Q) 7^ 0), we can apply further the duality 
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transformation on the k\ sum and obtain: 



x 



X 




3 /u*9»-l 

\ k 2 =0 

E 

V 



-«r( ^+2^fc 2 ) _ <7r jagazgg2±g*a fcl 

W2(wi w 2 — w ) 



^2^1 — ^-'^2+ a, ^2 



fci=0 



/^/ n \ / 1 -mA -™ — ^-^(e 2 -— ei) 
+ 0{1)\ — e "2e "i^2-^ 2V "2 

V ^2 

/■(".-£)«-' 
E 

y fc 2 =o 



(4.1.5) 



-iff " 2 a- fcf 2m 7 s 2-^2 fc 

g x g u^l^io^ — u J -)- (1/(1) 



/ 



Now taking into account the assumptions on u, the summation over the 

variable k x from to (uji — qi — 1 is the same as the summation on 

&4 = 0, . . . ,u\qi — 1. This fact occurs also for the fc 2 such that the above 
expression becomes: 



V' ' 01 (fi, 0) = ^ V'°' 0] (-fi"\ ir^H (4.1.6) 

(9(1) f 2 4°'°l (_q-i,q-^) + 2J°'?] 1 (-0- 1 ,Q- 1 ^) + lV 

where Pi, (resp. P 2 ) is the projection onto the plane {k\ = 0} ( resp. {k 2 = 
0}). We notice also that the duality transformation replaces essentially sums 
over k by other sums over k with an error of order 1. For example k2 is re- 

placed by (£ fel +(9(1)) (£ fc2 +0(1)) ~ E fcl , fe2 +0(1) E fcl +0(1) +0 2 (l). 
Moreover, there is a simplification since 

2 [o,o] f _o_i o-im _i Jo] f ^ utfz-uOA 

(4.1.7) 

For sufficiently small a;, integer sites such that < k\ < u±qi and < k 2 < 
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LO2Q2 are those with k = (ki, k 2 ) G O(A^). Thus: 



2J0.0] 



+ ©(lj.vg^c-n- 1 , 0-^ + 0(1)] (4.1.8) 

Now for a and 6 non zero (i.e. dj = {0, 1}; Oj = {0, 1}, % — 1, 2), we have in a 
similar manner the lemma 

Lemma 8 Vf2 G M 2 {jR) sym an invertible matrix and with u small enough: 

■ -iTr t en- 1 e _ _^ 

= 86 [>a™(-n-^-\e) +o(i).V^ tl (-n- 1 n- 1 ,g) 

+ 0(l).V^ Ad (-fi- 1 fi- 1 ,^) + 0(l)]. (4.1.9) 
wfcere = £ f e ^M)n(M) e ^-|)(M). 

4.2 Two-dimensional Gaussian sums with degenerate 
matrix Q 

In the expression of 2 <t^\Q,9) see Eq 14.1.11 we perform the duality trans- 
formation on k 2 sums. This operation yields: 



^292-1 t 2 

E 

fc 2 =0 



Sin7r gi (9, - f 2 (k 2 ~ 1)) i ^ i SmTTft^!-^ 



sin7r(0i-£(Afe-l)) sinvr^ 



(4.2.1) 



We observe that for Q non invertible, the general form of the sum 
e m quadratic form e m hnear form j g nQ ^ p reserve d un der this duality transforma- 
tion. Moreover the symmetry k\ <-> k 2 seems to be apparently broken since 
its expression depends on which variable to be transformed first by duality. 
This type of formula will be used later on, although with an appropriate form. 
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In the present case, it is sufficient to perform successive dualities uniquely on 
the variable k 2 until the summation on k 2 is equal to 0(1). Then perform at 
last the geometric summation over k\. We shall now remove the restriction 
to small uj and consider the general case. Let us give the following definition: 

Definition 9 An invertible symmetric matrix is called non- diagonal inte- 
ger, when there are no integer entries outside its diagonal. This set of such 
matrices is denoted by M*^(IR)i nVtS y m 

Let 

/ ^(jfe + fl) = e ^(fc+n)n(fc+n) e W(fc+n)e 

with 9' = 6 + Qn mod 1, we have 

= fmM^fm/^) (4.2.3) 

So for 9\ and #2 given, let (q[,q 2 ) suc h that the duality formula is valid, we 
cut the domain {0 < k\ < qi, < k 2 < q 2 } into small rectangles of size 
q[ x q' 2 . Let us consider the set: 

B(m 1: m 2 ) = {(k\, k 2 ) G IN 2 < < (nix + l)q[ and 

m 2 q' 2 <k 2 < (m 2 + l)q' 2 } (4.2.4) 

Necessarily m\ < e(qi/q' 1 ) and m 2 < e(q 2 /q' 2 ); here e(q/q') represents the 
integer part of q/q'. It is sufficient to verify that the duality formula is 
proved by pasting together the B(mi,m 2 ). Call m = (m 1 ,m 2 ). Then 

feens(o) 

— — * _ — * — * 

with 6 = 6 + Qm — N, where N is a vector with integer coordinates such that 
Vi, 0j = 9i+ (fim)j — A^j e]0, 1[. Now we perform the duality transformation 
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E W]W-W](^ det(n) 



-in t e'n.- 1 e' 



k<=B(m 



E /[-n-i^](ft (4-2.6) 



+ ( 1 )( E /[- n -sn-^](^)+ E /[- n -i, n -^](^) + 1 
y keVinB(o) kev 2 nB(o) 



with 



/[-n-n-^WW^™) (4-2.7) 



since *iV ■ m G N . Finally: 



^ /[n ' el(A;) " * v/det(Q) / [-^ 1 ^- 1 ^ (7V) 



k£B(m) 



+0(1) 



fee 

_!)fl(0) 



/ 



E ^'[-n- 1 ,n- 1 i9-f2- 1 iv](^) + E ^'[-n- 1 ,n- 1 e-n- 1 JV](^) + 1 

(4.2. 



fee 
•Pins(o) 



p 2 fis(5) 



\ 

/J 



or alternatively: 



E 

fc£B(m 



det(fi) 



/[-n- 1 ,n-i0](-^ r ) 



+0(1) 



E /[-n-i,n-ia|(^) 



(k-N)e 
QS(0) 



(fc-jv)e 



(fc-iv)e 
p 2 ns(5) 



/J 



(4.2.9) 



— * 

Here iV is an integer vector such that Vi = 1,2, |jVj — (fim)j| < 1. Thus 
replacing N by Qrh we make an error in the double sum J2 ki k2 of the order 
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0(1) ^2 kl and 0(1) Y,k 2 anci an error in res P ect to tne Zlfci ( res P- J2kJ of 
order 0(1). The index (k-Vtrh) G QB{6) is equivalent to k G H(m - 5(0)) 

— * 

or also k G fi.B(m)) by linearity. For the same reasons we have, for i — 1, 2, 
(fc — Jim) G V(£lB(0)) is equivalent to G V(VLB(fh)) and we have: 

91-1 qr 2 -l e(g 1 /q' 1 ) e{g 2 /q' 2 ) -i^SilS e ^l/q'i) e^/cQ 

E E /[n,$j(*0 = E E E Wl^^^/^fpy E E 



fc 1= k 2 =0 



r=0 s=0 



keB(r,s) 



r=0 s=0 



fenB(r,s) 



+ « E /[ (ft + E /[ 

\fcePinB(r,s) k&V 2 Q.B(r,s) 

E / - </!/••') (4-2.10) 

_feenA(gi,g 2 ) 



det(ft) 



+0(1) ( E /[-n-i^(ft+ E /[-n-^-i«|(ft + l 



Finally the same calculation shows that for all bounded domains V G IN 2 , 
we have the theorem 

Theorem 4.1 (Duality) Let VL G M£ (]0, 1 X> £/ie bounded domain in 
IN 2 we have: 



2\a,b] 



det(ft) 



VS-ir 1 ^- 1 *?) + o(i)/Vj™ (-n-sn- 1 ^ 



2 Ja,6] 



^i(nx>)v 



(4.2.11) 



This formula can be extended to the set of matrices f2 G M|(IR)i nt) . This 
result is given in the next section where the general case with d dimensions 
is treated. 
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5 d- dimensional Gaussian sums: Generaliza- 
tions 

In this section, proofs or parts of proofs that are similar to those in two 
dimensions will not be repeated. Let us consider the case where non-diagonal 
elements Wy, {i ^ j) are sufficiently small for the same vector 9, and where 

— * — * 

Q is invertible. We shall study the case a = b = 0. The general case can 
be deduced by analogy. Performing the dual transformation on only one 
variable k m with 1 < m < d on the sum defining: 



91-1 q d -l 

k 1= k d =0 
qi — l u m . m q m — l q d — 1 

= Yl ■ ■ ■ ■ ■ ■ Yl e l7rE1 ^ d " uk!+2 E i<^<™ ^ k ^+ 2 Zf<7< d M<] 

k 1= k m =0 k d =0 

q m -l 

X ^2 gWkm^ + ^^+EKKd"" 1 .)'™] (5.0.1) 
k m =0 



We get: 



91-1 (Jm-l 9d-l 

d af A (Q, 6) = ■ ■ ■ ]L • • • e i ^^ d ^ +2 ^^<^^ k ^ 

ki=0 km=0 k d =0 

X e 2 El<T<d k i * , / ! \ e -7J^( ™+52l<i< d Umiki) 2 



E 

x ( e -[-^+2(^+Ei< i < d ^fe i )^] + £,(!) 



f LOrnQm — 1 



km — 



(5.0.2) 



Let us define D m as the dual transformation operator on the index k m . We 
observe that for any function /: 

O m {l)f{. . . , k m , . . .) = 0(l)f(. . . , k m = 0, . . .) (5.0.3) 
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Thus regrouping terms using the previous remark, we have: 

— : 9 91 — 1 /Wmm?m-1 

I 



■Si=s m ( w ii-i-)*i ! 



fci=0 \ fc m =0 



^ e J7r (^m Ei^rJ"'" Wmin> K i) (5 4) 

fc d =0 

x e 2-ESr< d (^i-^^) feie 2 i7 rEgT< d 

Now this relation may be reexpressed using the dualization with respect to 
the variable fc m , T> m , defined by: 



(91-1 9 d -l 
fei=0 fc d =0 
91-1 /u> mm q m -l \ 9 d -l 

fci=0 V fe m =0 / fc d =0 

x e + 2 Ei<i<3<d^ fe < fe j+ 2 Ei<i<d fc «^] (5.0.5) 

— * 

Thus under the action of X> m , matrix elements of f2, are transformed and 
a multiplicative factor is introduced, leading to the following relations : 



J = Uij - ^£ii± vi < i < d, Vl< j < d with i ^ m and j ^ m 
j L_ 

A/ = g J^# m vi< i < d, with i ^ m 
6L = 



m UJ- 



mm 



$' = W— l — e l7r -mm$ with $ = 1 



and of course we have Q! = V m (Q), 9' = V m (6) and $' = £> m ($). Under the 
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action of duality on the index k m the matrix VL becomes a symmetric matrix. 
So let us apply this transformation successively as follows, 1 < m < d. 













— * 

9 








01 


_ 1 _ 









T> 2 



Or 



(5.0.6) 



where each bracket represents, with the set [Q, $] the effect of the trans- 
formation. Define now the symmetric matrix B m taken from Q as: 



LO U . . . L0 lr 



Ural 



UJr, 



(5.0.7) 



-iTr. t (e 1 ,...,e m )B m 1 (ei,...,e m ) 



Vdet B m 

with for 1 < i < m and 1 < j < m 



i,3 



det B m 



det 



^i— 1,1 



(5.0.8) 



LOij-i ^lj+1 • • • Wl,m 

■ ^i-lj+l • • • ^>i-l,m 

■ ^i+lj-l ^i+lj+l • • • u i+l,m 

^mj-l ^mj+i • • • ^m,m _ 



(5.0.9) 



that is essentially the (i,j) cofactor of B m . In a more compact way we can 
write 



&"m+l\(i,j)=l...,m — B m 



(5.0.10) 



for m < i and 1 < j < m 



m+l)i,j 



det B n 



det 



Um,l ■ ■ ■ <^m,j-l ^mi <^m,j+l ■ ■ ■ <^m,r. 



(5.0.11) 
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This is obtained by substitution of the j th column by the i th line. Hence if 
m < i and m < j, we have: 



1 



det B r 



det 



^11 • • • &l,m 



L0 ml 



(5.0.12) 



The right hand side is a (m + 1) x (m + 1) matrix obtained from B m having 
a piece of the i th line and j th column with m < i and m < j. Finally 



(5.0.13) 



and for all i > m 



- Qrn - 




' e 1 ' 


am 




_ dm _ 



®T — @i + y^(kWi)r>flj 

3=1 



(5.0.14) 



These results may be proven by recursion. The proof is uniquely based on 
the expansion of the determinant into minors. For example for (a; m +i)n, one 
has to compute: 



(CO 



(^m)ll — 



mjlm 



(^m)r 



(5.0.15) 



with 



(Wm) 



11 



(B m )u 

det Brr, 



where B m is the matrix formed by cofactors of B ri 

det B m+ i 



det B r , 



(5.0.16) 



(5.0.17) 



With det B m+ i — LU m +l,m+l det B m + Y1T=1 Y^j=l W i,m+l w m+lj(^m)ij 
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Ui, m+ i det "y ' (5.0.18) 

Hence 

/ , \ Ki)lm ^m+l,m+l(-Bm)ll /'cnm^ 

(^m)ii - t — s — = j— -= (5.0.19) 

\^m)mm Get %+l 

j mm 



det 13,,,-- , . . 

1=1 3=1 



Knowing that the determinant is invariant (up to a sign) under permutations 
of lines or columns, it is sufficient to compute: 



m—l 



(B m )n(B m ) mm — (B m )lm(B m ) m l — I ( — Y\ w i+l,(r(i)+l 

y<re5 m -i 1=1 

m—l \ / m—l 

e (-D' M n --'o - e (-D <M n^(.) + . 



^'sSm-i i=i / \o"e'S m _i i=i 

m—l 



e (-i)^ n ( 5 -°- 2 °) 

v cr'e5 m _i 1=1 J 

Where S m -\ is the symmetric group of permutations of order (m — l). After 
expansion and regrouping, we get the double sum: 

m—l m—l 

E E (-iy {a)+e{a,) Yl II (^+l,a(*)+l^>0) - ^^(ij+iWj+i,^-)) 

o-e5 m -i cr'e5 m _i i=i j=i 

(5.0.21) 
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Now the product of permutations a <g> a 1 G S m -\ x «S m -i may be represented 
as crCSo"' € <S m x iS m _2 and the previous double sum is recast under the form: 

m m—2 

e e (-i) fW+fM nn(^)^ + i,o- )+ i) (5.0.22) 

<reS m o-'e5 m -2 »=1 j=l 



i,<r(i) ' 



m-2 



= det _B m det 



^12 



<^l,m-l 



^m-1,2 • • • ^to-1,to-1- 



For opposite permutations of lines and columns we have the following rela- 
tion: 



{Bm)ij{B m )kk' — (B m )ik(Bm)jk> — 

Uli . . . U>ij . . . U>i ; jfc' 



det B m det 



CO 



hJ 



CO 



(5.0.23) 



Wkl ■ ■ ■ Wkj ■ ■ ■ U)k,k' ■ ■ ■ ^k,r, 



(5.0.24) 



The last matrix is of size (m — 2) x (m — 2) obtained from B m by removal of 
columns j and k! and of the lines i and k . Consequently: 
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con 
tin 
tiki 



ML 



/ o \ mm 

det£ m+ i V^ft^ det^+i 6 



a;. 



Wl,r, 



\ 



(5.0.25) 



LU m \ . . . ti m j . . . ti rn ^i . . . ^m,m _ ) 

But this expression corresponds to the expansion of the determinant of 



(^m)ll 



m+l)U 



det B m+1 



det 



^22 



Wro+1,2 



<^2,m+l 



<An+l,ra+l 



(5.0.26) 



with respect to line (m + 1) and column (m + 1). Other results may be 
proved along the same line, i.e. using the expansion of the determinant into 
minors. Hence if Q is invertible, one can perform the duality transformation 
with respect to all indices kj, and the constants have the expressions: 



d 



V? 



t en- 1 e 



Vdetfl 



(5.0.27) 



Each sum J2 k , is replaced by J2k t +Ci(l)> as seen before. Proceeding as for 
d = 2, one can extend the two-dimensional result to the set of invertible 
matrices with non-integer diagonal elements. Due to linearity, this formula 
remains valid for sums on a bounded domain T> e 7L d . The hypothesis of 
a bounded domain is necessary since otherwise the Gaussian sums are not 
converging. This is not the case for theta functions. Consequently, the result 
is also valid for matrices il G M d (]0, 1 

There remains the extension to IR \ TL to be examined. For this, let us 
introduce L^ e M d ({0, 1}) such that (Lij) k ^/ = for all k,k' except for 
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(k — i,k' — j) or (k = j, k' — i). Compute now: 

d o^\Q, + Lij, 9) = £ e KHM)<^)(M>+*<*-i>(*-§>) (5.0.28) 

kev 

= iaiaj ^ e ivr(*(n-|)n(n-f)+2 t (n-f)(e-|)-(a l n J +a,n l )) 

kev 

since: 



(n - -)^(n - -) = 2(n^ - — )(nj - —) = 2n inj - (a inj + a^) + } 

(5.0.29) 

We set b' = b + Oje} + mod (1, 1) and obtain: 

'o**(a + L„, «) = (-,:)--' £ ^('C-IW-IH'C-IM*-*) 

fee© 

= (-i) a ^ d a ^ ff] (n,9) (5.0.30) 
Now let Li e M rf ({0, l}) sym an d {L>i)k,k' = except for k — i, k' — 1. Then: 

fee 

(5.0.31) 

— # — # 

but with 6' = b + (1 — a^ej, we have the result 

d a%'\n + Li,&) = (Vi) aid <J [ ^ ] (n,9), (5.0.32) 
and may state the theorem: 

Theorem 5.1 With the convention = i and \/Q e M^([0, l[)i nv and 

\/9 e]0, l[ d ; £/ie Gaussian sum in d dimensions admits the duality formula: 
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v v ; detfi ^ vl 1 ...v % d d sm K ' 

0<il,...,i l l<.d 

(5.0.33) 

We use also the convention V° = Id, the identity map. 

• For Lij G M d ({0, l}) sym introduced before and b' = b + a^e} + aj€l 
mod (1, 1) we also have: 

<Vg' S] (ft + L^, 8) = (~t) a ^ d a§' bl] (n, 9) (5.0.34) 

b + (1 - ai)ej, 

(5.0.35) 

The functions d C^ b \Vt,0) and d S^ h \Vt,0) verify similar duality relations. 

If Q is not invertible, one may perform a partial duality transformation. 
Up to a permutation of indices, the duality transformation may be applied 
to those of the indices for which the minors are non-zero and the number of 
indices corresponds simply to the rank of the matrix Q. For the remaining 
indices the summation can be performed, since they are geometrical sums. 



• And for Li G M^({0, l}) sym introduced before and b' = 
we obtain: 

d a [ i' $ \n + l u 9) = {Vir d a [ f\n, 9) 



6 Applications to the computation of £ n [f3, X] Pi 



Now back to Eq. 13. 0.71 which we shall write using the notation of the previous 
section. The matrix Q belongs to M^^IR)^™ and has the form 



n 



ft 



p 2 (E^-^) -P2P3 
-P2P3 p 3 (Y,Pi-P3) 



—?2P2q 
-P 3 P2 q 



-P2?2q 



-P 3 P2 q 



(6.0.1) 
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we note E for Ei=2 an d 



e = x 



P2 
P2, 



e Vect 2q -i(fR) 



(6.0.2) 



It is easy to see that det £7 = (rank(fi) < 2(q — 1)). By just adding all the 
lines, we get a line of 0. Let us compute the minor: 

p 2 (J2P i -p 2 ) -p 2 p 3 ... -P 2 P 2q -i 

-P2P3 Ps^Pi-Ps) -P^2 q -l 



det 



g 



— P2P2q-\ 



3 2(q-l) P 2 ■ ■ ■ P 2q-1 . 



Y,Pi-P2 ~ P 2 ••• -^2 

-P 3 : -P 3 



2g-l 



°2q-l ■ ■ ■ E "^i — P2q-\ J 

(6.0.3) 



Now replacing the I s line by the sum of all the lines, we obtain: 



f3 2 (i-Vp 2 ...P : 



2q 



( E p t )2( 9 -1) 



det 



1 1 ... 1 

-p 3 y i', r, ; -p 3 



—P2q-\ — -P2Q-1 ••• E^ - -P2Q-1 

(6.0.4) 

Replacing the j th column by the difference of the j th and first column. Doing 
it for all j > 2, we are led to the expression: 



_ g 2(q-l) P2---P2q , 



1 











-p 3 


ZPi 


















--P2g-1 











= I3 2 ^ 



i) 



P9...P 



2q 



(E^) 



7^0 



(6.0.5) 
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Consequently we have the lemma 
Lemma 10 

rank(Q) = 2{q - 1) 



(6.0.6) 



This result implies the duality transformation of £ n only with respect to 
indices k 2 , . . . , k 2q -\. Thus we have: 



2(9-1) 



(9-1) 



P-1...P, 



i e -i^ t (62,-,e2 q -i)o 2( 1 q _ 1) (e 2 ,-,e 2q -i) ^ q ^ 



2q 



with 



d lf o _ 1 

i,j=2,...,2q-l — It — 



p 2 (J2P i -p 2 ) ... -p 2 p 2 ,_i 

— PlPlq-l ■■■ p2q-l(52 Pi ~ ^2g-l) 



(6.0.8) 



Using the results of the above section we also obtain Wi = 2, . . . , 2q — 1 

1 



w 2g-l|i,2g 



det 



p 2 (^p i -p 2 ) ... -p 2 p 2q 
-P2P3 



(6.0.9) 



-P2P 

— -P2-P 2g -1 



— P2qPi 
— P2q-lP2q 



— P 2 P 2q -l 

— PlPlq-\ 



— P%P2q-\ 
~P2q-l(Yl Pi ~ P^q-l) 



In the last matrix we substitute the i th column by the sum of all the columns 
of the matrix such that i th column becomes: 

[~PiP 2 , P l {P 1 •...•/';• ...P 2q ) -PiP 2 q-l] (6.0.10) 

Such that u 2q ~i\i t2q = l.We also have: 







02 




= Q- 1 








02q-l 


n2q 
U 2q- 


1 = 02q + 


...+02 



(6.0.11) 
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At the end we get the form of il 2? - 2 : 



&2q-2 — 



n- 1 
i . . . i 



i 
o 



(6.0.12) 



where is the (2q — 2) x {2q — 2) matrix defined above. So considering all 
the results modulo 1, the matrix f2 2 <?-2 does not admit coupling between fc 2g 
and ki for i — 2, . . . , 2q — 1. Consequently with 



Afj, — A 



n — 1 n — 1 



2q — 2 times 



and #i = rrPi, Mi = 2, . . . , 2q - 1 



£> 2(? -i o • • • o £> 2 (£ n (/3,x;P 2 ,...,P 2(? )) 







9-1 



P 

p 2 ...p 2g 



(6.0.13) 



n— 1 n— 1 

e ^'(ea,...,^,-!)^- 1 ^,...,^,-!) ^ ^ g27r:, ■/,•,,,-/ I /' 



0(1) 



i2+- +i2 9 -i 2q-2 [0,0] 



E 

fci=0fc 2 q=0 2<i 2 ,...,t2q-l<2q-l 

Q- 1 Q- 1 



v P ' P 



-0 



where stands for the vector: 



= 



02 
02q~l 



(6.0.14) 



Hence: 



/ \r^2q n-1 n-l 

UP, x;P 2 ,...,P2q) = 0(l) J J2 E e 2 ^ k ^ S 



^ P, 2 o...oP i25 _ l[(5As] \^ > J 



(6.0.15) 



2<i 2 ,...,i 2q -i<2i}-l 
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Here we are reduced to compute the a-functions, with f2 a matrix with integer 
entries of the form: 



a 



[0,0] 
Aft 



xn, if) = J2 e 



in(\ t knk+ t 



(6.0.16) 



where A 6 1R 



— * 

Lemma 11 For any matrix Q with integer entries, there exists a basis (k 1 ) 
such that the quadratic form l kVtk becomes: 



kVtk = \ik[ 2 + . . . + \2q-2Kq-2 



t k'Ak' 



(6.0.17) 



where A is a diagonal matrix with rational entries. If we note k' = Uk, U 
has also integer entries. 

Proof. Let us write down explicitly the expression of the quadratic form: 



^ 1 ^ijkjkj — ^ ] uJuki +2 ^ ] LUijkikj 

l<i<j<d 



Ki<d 



k m + ^ ] 



U), 



-ki\ - LO r , 



=1, i^m 



E 

\i=l, i^m 



UJ r , 



Uijkikj 



+ 2 £ 
= E ^ + 2 E 



uj'^k'ik'j (6.0.18) 



i=l, i^m 



l<i<j<d, (i,j)y^m 



So by redefining the u^- and proceeding by induction we obtained the ex- 
pected result. ■ 

The transition matrix U used for the change of basis can be taken with in- 
teger entries instead of rational entries. Thus Vi = 1, . . . , d, k\ — X^=i Bijkj 
and U is of course invertible. Moreover, up to a permutation of indices one 
can insure that the next minor is non-zero. We set then 



U' 



U d - hl 



Ud-14-1 



(6.0.19) 
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and det(U') ^ 0. Let V be a bounded domain of IN d and let P G U{V). Then 
P, by construction, is a vector with integer components and 3k <E V such 
that P = Uk. At this stage, one should show that U(T>) can be decomposed 
into integer sub-lattices Ri of lattice spacing Si G IN. Thus P G Hi, then 
P + bSj G Pj with efj, j — 1, . . . , d is a unit vector. Pj is defined inside the 
boundary of U(V). To prove this last property one must solve in IN d , the 
system Mi — 1, . . . , d — 1: 



3=1 



We have thus 





h 




u ltd 


u' = 




= ~k d 






kd-i 







and as det U' ^ then: 



fci 

k d -i 



kd _ & 



det P' 



P M -i 



(6.0.20) 



(6.0.21) 



(6.0.22) 



where U' is a matrix of cofactor of U' (with integer entries). Let us take 
k d = — det U' so we have determined in !N d the solution the above equations 



h 








= U' 











(6.0.23) 



It is clear that any multiple of this vector is a solution. Let us compute now 
S, defined by: 



d d-l d-l 

5 = Udjk, = U dd det U' + Yl UdjU'ijUid 

3=1 i= l 3=1 



(6.0.24) 
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Thus 5 = det U by expanding with respect to the last line and column, we 
obtain: 

Proposition 12 Let U(V) be a bounded domain of IN 1 and U an integer 
invertible matrix . Then there exists at most 5 d hypercubess Ri G IN 1 of lattice 
spacing Si (Si being a divisor of 5, with S = detU), such that VP G U{V), 
31 < i < S d for which P G Ri. The hypercubes (Ri)i=i,...,s d are a ^ limited by 
boundaries of U (D) . 

Proposition 13 Let (^ 7 n) ng ^v be a family of hypercubes of IN 1 , defined by 

— * 

k = (ki, . . . ,kd) G T n with VI < i < d, ki G IN and 1 < ki < n. Let ("Cn) ne ^y 
be an increasing sequence of connected domains of IN 1 such that there exists 
t G IR for which Vn G IN: T n C T> n C T tn - Then, uniformly with respect to 
9: 

d 4^6) ~ d afSW) (6-0.25) 



n— >oo 



Proof. 

IN d as a partially ordered set, we can easily prove that V(a, b) G Vect^O, 1}) 

, Vfi G M d (JR) and V# G Vect d (5R) the mapping d a^j b \tt, ff) defined as follows 
is increasing: 

d af^9): V{^ d ) — C 

j= <Vg' bl (Q,0) 

So for n sufficiently large we have: d a l ^ ] {Q,ff) < d a l ff{Q,6) < d a [ ^{Q,6). 
Now taking in account that: d af b }{Q,6) = 0( d a [ ^\n,9)), we can conclude 
that d <j [ vf{tt,6) = 0( d af n b] (n,e)) and is homogeneous in 6. ■ 

Remark 14 

Here we just recall the result of [DM] that: Vx G [1/2, 1] 3x G IR such that 



According to [HL pp. 202], the hypothesis of uniformity is automatically 
verified in Gaussian sums with d — 1. This hypothesis remains valid for 
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ci-dimensional Gaussian sums. Let i?i(n), . . . , Rd(n) the sub-lattices making 
up B(QAfi). Then we have: 



l n-1 



ir 1 n 



0<i 2 ,...,i 2q -2<2q-l H y 

V V ofi^+--M>.-i 2 9 - 2 [o,o] f _A 1 _A l 3s 

j=l i 2 e[i,2 9 -i],..., 

»2,-2£ [1.29-1] 

(6.0.27) 

with 0' = B~ 1 6, where B is the matrix of basis change described in lemma fTTI 
such that A is a diagonal matrix with integer (or rational) entries. Since A has 
integer matrix elements, it follows from the above formula that there exists 
a finite number of sub-lattices for which the multi-variable sum is reduced to 
a product of (2q — 2) factors of one- dimensional Gaussian sums where each 

factor is of the form 1 o r [ ( ^ |(— A, — jpO) where A is an eigenvalue of A. The 
following section will consider the problem of the asymtotic behavior of such 
Gaussian sums. 



6.1 Estimations of Gaussian sums using the duality 
formula 

We shall now present the theorem which motivates the previous construction. 
Let us recall the method of Hardy et Littelwood using the duality formula for 
the computation of the Gaussian sums. They give the expression, Vx g]0, 1[, 
V0 e]0,l[ and V(a, b) E {0, 1}, and any n E IN: 

1 a^ b \x,B) = 0{l) y /nxx 1 ...x v + °^ -. (6.1.1) 

y XX\ . . . x u 

Recall briefly the steps: For convenience let us denote xq = x. Vz = 
0, . . . , v, we have < Xi < 1. Xi are the rest of the development in continuous 
fractions of x: 

x= x n = [oi,o 2 , ■ ■ ■ ,] (6.1.2) 
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and Vz > the (a?i) ig ]N are given inductively by the relation: 

Xi_i = — -J- — (6.1.3) 

Hence, for any n, there is an integer v determined in such a way that : 

nx . . . x v < 1 < nx . . . x v _i (6.1.4) 

Through this condition, the first term on the r.h.s. in Eq.6.1.1 becomes 
irrelevant with respect to the second one. Let us put v = u(n) in the above 
expression. So there exists a positive constant H such that: 

x . . . x v > Ha~( n) x . . . x v - X (6.1.5) 
Hence we obtain by using Eq J6.1.4l and Eq 16.1.51 

1 0(y/a^n) (6.1.6) 



y/X ...X„ 



If we use the trivial estimation (xiX i+ i < |) we deduce that necessarily 
Ve > 0: 

2 + e 

< ^ - Inn (6.1.7) 

thus: 

• if a n = 0(1) then l at b] (x,9) = 0(^/K) 

• if a n = 0(nP) then l at b \x, 6) = 0(Vn(lnn)§) 

• if a n = 0{e np ), with < p < ^ then l at b] {x, 6) = 0(^/nn& +e ) 

Nevertheless let us recall that the equidistribution of e l7rxfc2 , for any irrational 
x, [K] and [KN] give the upper bound of: 

1 aM{z,0) = o(n) (6.1.8) 

Now we must use other results of number theory concerning the computation 
of continuous fraction expansion of a fraction of the numbers occurring as 

arguments in our Gaussian sums 1 o r {'^°j(— jp, ~Jj3@) > °f the form: 

^frr-^^-W 1 ' )) (6 - L9) 

V tx + w I 
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where (u,v,t,w) G IN 4 . 

We will refer to the article of Raney [R] and to the survey of Van der 
Poorten [VdP]. Let us introduce some notations: 



M> fe : ]0,1[ — 

x i — ► [ai,a 2 , ■ ■ ■ ,a k ] 



such that 



(6.1.10) 



ai + 



02 - 



(6.1.11) 



Denote \I/oo the map which associates to any real number its continuous 
fraction. is a homeomorphism from ]0, 1[ to IN*^. Referring to the 

theorem of convergence for continuous fraction, for any x e]0, 1[ we have: 



X 



(6.1.12) 



Consider now the matrix representation {R, L} of the expansion in continu- 
ous fractions of x: 



R 



Here Wk G IN* one gets 



R k 



"11" 

1 


L = 


"10" 
1 1 


'lk' 

1 


L k = 


" 1 

k 1 



So ^ k 1 ([ai, . . . , a k ]) G Q is determined by the formula: 

W fc Qai, . . . ,a k \) - ^ 



E'=i(^ ai ^ a2 ---^ afc )2, l 



(6.1.13) 



(6.1.14) 



(6.1.15) 



Now we use the main result of the articles quoted above: 

The continuous fraction of is related to that of x if ad — ac ^ 0. To 
this end, let us consider the set of 2 x 2 matrices of positive integer entries: 



M 



a b 

c 



b d ^j G M 2 (IN)| n = detM ^ 0, a > c,b < dj (6.1.16) 
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The number of such matrices is finite [R]. Denote by S n = {Ai, . . . ,A N }. 
The commutation relations obey some restrictions: Mi = 1, . . ., N, 3i ,i" = 
1, . . . , N such that: 

AiR k nL k L = L^R^A, 

AiL k *R k R = R^L^A,, (6.1.17) 

with k\) positive exponents such that 1 < k l R + k l L < n and so on. The 
algorithm which allows to compute, for a given n, the matrices in the left 
hand side is given below. We note that these commutation relations are 
invariant under transposition: 

AiR^L^ = L k ^R k RA,., 
t A i 'R k R Lfr = L^R^Ai (6.1.18) 
Note that for any n G IN* 



n 

1 

n 
1 



R = R n 



n 

1 

n 

1 



(6.1.19) 



This provides an algorithm to compute the expansion into continued fraction 
of the values of the function x i— > 

cx+d 



6.2 Illustration of the algorithm 

We now shall illustrate this method for n — 2. The set £ 2 is made up two 
matrices which represent the multiplication and the division by 2: 

£2 = {A, A'} (6.2.1) 

- 4 =[o?] A ' = [ll] (6 ' 2 - 2) 

The rules of commutation are compiled in the following tabulation: 





A 


A' 


A 


R : R 2 
L 2 :L 


LR : RL 


A' 


RL : LR 


R 2 : R 
L : L 2 



(6.2.3) 
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The table is to be red following the exemple of the first entry:v4i? = R 2 A. 
These relations are equivalent to the following distinct cases: 

x = [0, 2a, b, . . .) 2x = [0, a, 2b, . . .) . . 

x = [0,2a + 1,6,...] 2x = [0,0,1,1, (6- l)/2, . . .] 1 ' 

The entire determination of 2x from x is done by induction with respect to 
those rules. For example, if x = [7, 2,5,.. .] then 

2x = [3, 1, 1, 0, 1, 1,2,...] = [3, 1, 2, 1,2,.. .] (6.2.5) 

Let us introduce some notations: [x , Xi, . . .] with < Xj < 1, Vz G IN, the 
sequence of the rest deduced by the decomposition into continuous fractions 
of x = x . We shall denote the resulting sequence of the rest of 2x by 
2x — x' — [x'qjX^, . . .]. Thus we have , by definition, the relation: 

Xi = -b i+1 + — (6.2.6) 

where x = [0; b\, 62, . . .]. As the formula Eq 16.2.41 shows, the processus of 
commutation (resulting of the multiplication by 2 generates a lag length 
between the two sequences (2 = [0; 61, b 2 , ■ ■ ■) and 2x = [0; b^, b' 2 , . . .]). We 
shall define a function (depending on the decomposition of x) which measures 
the length of 2x with respect to x. 

K -™^ KM ™ 

such that : 

• if 2 divides 61 ( in what follows denoted as: 2 | 61 ), then k' x {ki) = k x 
and k' x (k 2 ) = k 2 

• if 2 does not divide 61 ( in what follows denoted as 2 { 61) then k' x {k\) = 
ki + 3 and 

- if 2 f b 2 then k' x (k 2 ) = A4(/ci) + 1 

- if 2 I b 2 then k' x (k 2 ) = k' x (h) + 3 

and so one... We recall the trivial reduction which occurs when one of the 
components of the decomposition is zero, we shall refer to the example Eq 
Note for this x , ^(3) = 5. 
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Lemma 15 Vp & IN, let [xo, x±, . . .], [x' Q , x[, . . .] be the rests of the expansion 
of x and x' = px respectively, then Wk e IN 3C G M which divides k (possibly 
C = 1). such that 

x i--- x 'k' x (k) = Cx i ■■■ x k (6.2.8) 
where [x , Xi, . . .], [x' , x[, . . .] are respectively the rest sequence of x and px. 
Proof. 

We shall reduce the proof to the simplest case p — 2. The other case is 
quite the same but it turns out to be more difficult because of an involved 
algebra. We prove the lemma distinguishing all the possibilities: 

• if 2 | b\ then b[ = y. Compute now x' x' 1 = x' ( — ^i + -^rj = 1 — 
bix = xqX\. Consequently x[ = y. Since 2 | b\ then b' 2 = 2b 2 and 
x\x' 2 = 1 — b' 2 x[ = 1 — b 2 x\ = x\x 2 and x' 2 = 2x\. 

At this level we have two possibilities either 2 | 63 and the procedure 
repeats itself identically or 2 \ 6 3 , in this case we shall refer to the 
following 

• if 2 \ bi then b[^Y^ an d b' 2 = b' 3 = 1. Hence let us calculate XqX^x^x^: 
x'qX^x'zx'z = -(26i + l)x' Q + 2 = -2b Y x + 2 = 2(1 - Mo) = 2xix so 

We establish another relation x'qx' x x' 2 = (1 + b\)x' Q + 1 
then x' x x' 2 = 1 + b\ — jj- = ^p- ie x 1 + 2x' 1 x' 2 = 1 and we deduce that 
(1 — Xi)x' 3 = 2x\ 

• If 2f& 2 then& 4 = ^=1 an d6' 5 = 2b 3 . Calculating x' 3 x' 4 x' 5 = x' 3 (l+b' 5 b' 4 )- 
65, so 0:4X5 = 1 + b 2 b 3 — ^ = x 2 x 3 . We also prove that — b' 5 x' 4 + 1 = 
X2X3 = —b?,x 2 + 1 so 14 = ^ and x' 5 = 2x 3 . 

• If 2 I b 2 then b' 4 = ^ end b' 5 = b' 6 = 1. x' 3 x' 4 x' 5 x' 6 = 2 - (b 2 - l)x' 3 . So 
using 2x\ = £3(1 — xi) we get x' 4 x' 5 x f 6 = x 2 . we also get 2x' 4 x' b + x 2 = 1. 
So let us summarize these results in the following tables: 
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X = 


2b x 


6 2 


2b 3 






Xn 


X\ 


X2 


X<\ 






2x = 


h 


2b 2 


h 






x' Q = 2x 


rr.' — XI 

x l 2 


x' 2 = 2X2 


rr' — %3 

x 3 ~~ 2 






X = 


26i 


6 2 


26 3 + 1 






Xq 




x 2 








2x = 


h 


2b 2 


b 3 


1 


1 


x'q = 2xq 


rr' — XI 

X l — 2 


x' 2 = 2X2 


4 

<-y» <T™ r Y* T* 

x 3 x 4 + 2x 3 = 1 


x^ 


X 5 



(6.2.9) 



and we also have: 



X = 


26i + 1 






262 + 1 


63 




Xq 


Xi 








^3 




2x = 


bi 


1 


1 


62 


26 3 




x'q = 2xq 


X 

iy* i~rJ rpi rpf 

Jo y — «x.' ^ tXj 2 ^ 

x'^2 + 2xi = 1 


x 2 


x 3 


2 


2x 3 




X = 


26i + l 






26 2 + 2 






X Q 


Xi 












2x = 


61 


1 


1 


6 2 


1 


1 


x'q = 2xq 


x\ 

x'^2 + 2xi = 1 


x 2 


x 3 


— X^X^Xq 

X4X5 + 2x 2 = 1 


X 5 


x 6 



(6.2.10) 

Thus we have : \/k e TN x± . . . x k = x\ . . . x' k ,^ or 2xi . . . x k = x[ . . . x' k ,^ k y 
This complete the proof of the lemma. ■ 

More generically, let A E S n goes through the string R ai L a2 R a3 . . ., com- 
muting with respect to the two matrices {R, L}: [AR ai L a2 R a3 . . . R ak ] 1— > 
[R a iL a 2R a 3 . ..R a k'A'} with A' G £ n . This procedure gives an algorithm for 
the computation of: 



(6.2.11) 
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which only depends on [ai, 02, . . . , a^]. Now using the above lemma in the 
formula Eq I6.1.l[ we obtain the theorem: 

Theorem 6.1 V(u,v,t,w) G IN 4 , Vx G M \ Q, V0 G [0,1[ and V{a,b) G 
{0, 1}, ffcen 

V n' 6] (^^)=°( V " ,Q] ^°)) ( 6 - 2 - 12 ) 

By applying this theorem to the above expression one can deduce that: 

7t — 1 n— 1 



e 



2l7ra(fc2g-fcl)(E'Pi) 



fcl=0fc 2 g=0 



x C(V[ a ' bl (/5,0)) 29 - 2 (6.2.13) 

Remark 16 

The estimation depends on (u,v,t,w). It remains an essential question: 
does there exist a constant K independant of (u, v, t, w) G IN 4 such that 

vr ] (^^)<^.(vr] (X)0 )) (6.2.14) 

This result is not trivial to prove. In fact, if this assertion were true, the 
characteristic function would be analytical in a neighborhood of zero. Con- 
sequently, this work only gives a partial answer. 



Remark 17 



We observe that this result is an optimal estimation [HL p225, theorem 2.221] 
when applied to the matrix case: for any A is an integer invertible matrix, 
and for any (3 G IR \ Q, we have: 

d a [ ^ ] ([3A,6) = OCa [ y ] ((3,0)) d (6.2.15) 
It cannot be replaced by a better estimate. 
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Now we can compute the mean value of S^ q (^) (see Eq. IH.0.5J) : 



E 



Pi 



(6.2.16) 



where J} stands for rii=2- Performing the integral in x yields the condition 
ki = k2 q , we end up with: 



E 



^(Vr'^O)) 2 ^ 1 ) 

(n^) | v / E^ 



(6.2.17) 



O^ai? (x, O)) 2 ^ -1 '' depends on n but also on (P 2 , . . . , p2( q -x))- Moreover 
Cqis finite because the moments are perfectly defined (the integration is done 
over a compact set [0, 2n[x [0, 2tt[). So: 



E 



n.0 (Vj^ (/3,0) 2(9 - 1} 



(6.2.18) 



Remark 18 

For g = 1, we recover the exact result E [Sl q (^f )] = which is independent 
of the choice of (3. 

Theorem 6.2 T/ie unique sequence f n ( up to an equivalence as n — > oo) in 
order that: 

2q 

= 0(1) (6.2.19) 

with respect to n for all q, is f n = -y/n. This estimate is true for all (3 such 
that l a [ n'° ] {(3, 0) = 0{y/n) 



E 



Mf) 

fn 
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Proof. 



We must choose a (3 having a bounded continuous fractions representation 
to have /„ indenpendent of q. The proof results straightforwardly from the 
remark [TH that is: 



Remark 19 

Note that if (3 <E Q then [HL] VL°'° ] (/?, 0) = 0(n) and we recover the 
result established earlier, i.e: 



which means that there is no normalization of S n for which the moments 
neither diverge nor be non-zero. 

Examples of numbers (3 which admit an expansion as an infinite and 
bounded continued fractions are the quadratic irrationals. They have ul- 
timately periodic continued fraction expansion. This class contains all the 
square roots of products of pairwise distinct prime integers, and of course we 
have, for those numbers, Vg G IN 



Hence there exists here a normalization ^Jn of S n for which all moments 
remain bounded with respect to n and are not zero. 

Remark 6.3 

Any attempt to obtain information to the convergence in distribution of 
S n / fn by using numerical computations cannot be correct since behaviors are 
radically different according to the nature of the number (3 and the rational 
numbers are dense in IR. Finally all of these calculations may be applied in 
the same way to any periodic function instead of the signum function x{y)- 
As to the problem of the convergence in distribution of S n / y/n, we shall 
make some comments in the conclusion: 




(6.2.20) 




(6.2.21) 




(6.2.22) 
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7 Conclusions 

In this approach to determine the limiting distribution law of we have 

sought to compute the asymptotic behavior of the moments of If the 
parameter (5 — ^ admits an expansion in bounded continued fraction, the 

behavior obtained for the expectation values of moments E ( ) = 0(1) < 

Ak may lead to the convergence of the series Xlfclo ^lr^ fc ' aroun d t — 0. 
In that case, the sequence $ n (t) converges towards $(t), which is analytic 
near the origin. This implies the existence of a limiting distribution with 
finite moments. The estimation of the speed of the increase of the Ak seems 
difficult and is still an open problem. It is thus difficult to have an idea of 
the limiting distribution without such estimation, Although the procedure 
does not imply the convergence in distribution for it shows that this 
normalization leading to bounded second moment, fails to lead to bounded 
moments of higher orders when (3 is an irrational having no expansion in 
bounded continued fractions. 
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